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Abstract. A novel approach to analyzing time series generated by complex systems, 
such as markets, is presented. The basic idea of the approach is the Law of Self-Similar 
Evolution, according to which any complex system develops self-similarly. There always 
exist some internal laws governing the evolution of a system, say of a market, so that 
each of such systems possesses its own character regulating its behaviour. The problem 
is how to discover these hidden internal laws defining the system character. This 
problem can be solved by employing the Self-Similar Approximation Theory, which 
supplies the mathematical foundation for the Law of Self-Similar Evolution. In this 
report, the theoretical basis of the new approach to analyzing time series is formulated, 
with an accurate explanation of its principal points. 
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1. Law of self-similar evolution 

There are two sides of the approach I am going to present here: philosophical and 
mathematical. The first side gives the general idea of why should we think in these or 
those terms, and the second side is supposed to answer how concretely could we realize 
the idea. It is natural to start with the formulation of the general idea. 

The idea to be formulated is general since it may concern any complex system, for 
instance, a market, a firm, a society, a nation, species, a person, a brain, a heart, and 
so on. Shortly speaking, a complex system is such that cannot, at least at the present 
time, be completely described by any finite number of given equations. A complex 
system exists as an empirical organism, with its specific character and behaviour. 

To understand the character of a system means to notice some rules governing its 
behaviour and, therefore, to be able to predict the latter. If the system is complex, 
such rules cannot be described by fixed equations. But one may always characterize the 
behaviour in terms of self-similarity, implying that there exist some specific features 
that, to some extent, are repeated during the evolution of the considered system. Thus, 
characterizing a person, we name the corresponding features that give us impression 
of what can be expected from him or her. 

The evolution of each complex system is always governed by some laws. This is 
why there is a similarity in the behavioral facts of this system. Such a behavioural self- 
similarity is easily noticeable in the life of particular persons as well as in the history 
of societies, nations, and biological species. This observation suggests that there exists 
the Law of Self-Similar Evolution, according to which the main features of a complex 
system are preserved in the course of its evolution. And if one has grasped these 
features, he or she should have the feeling of what could be expected in future. 

When the behavioral facts, related to the evolution of a complex system, can be 
associated with some quantities, one gets time series. These can be market prices or 
indices, measurable signals from a brain or from a heart, some data characterizing 
societies or species etc. The examples of time series are numerous, their consideration 
being usually based on statistical analysis and the construction of dynamic regression 
models [1-5] . These models provide a rather reasonable description of sufficiently stable 
dynamics with repeated events, like seasonal variations, but they fail in treating such 
unstable systems as markets. Financial markets are very difficult to predict. Perhaps 
the most central question in finance is under what circumstances is prediction possible 
at all? 

The widespread opinion is that market prices fluctuate absolutely randomly and, 
thus, are unforecastable. This point of view is what is called the Efficient Markets 
Hypothesis, which can be traced back to Samuelson [6]. According to this hypothesis, 
in an informationally efficient market, properly anticipated prices fluctuate randomly. 
This includes the assumption that prices fully incorporate the expectations and infor- 
mation of all market participants, that all market traders are identical in the sense that 
all of them possess the same and the whole information on a market, all traders having 
the same physical and mental abilities. Such a hypothesis that all investors are fully 
rational agents that instantaneously and correctly process all available information is 
clearly unrealistic. Moreover, people, thanks God, are never identical neither in their 
abilities nor in their wishes. They are not elementary quantum particles, but each of 
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them is a very complex system by its own. There is also increasing empiric evidence 
that even the most competitive markets are not strictly efficient [7-9]. 

A word of caution is called for with respect to the meaning of the term "random", 
which is often confused with "chaotic". These two notions are rather different. A real- 
istic process can be chaotic, at the same time, possessing properties of both regularity 
and randomness in different proportions, so that it is, in principle, possible to distin- 
guish chaos from randomness, and even control or erase chaos and make predictions 



Dynamics of all realistic complex systems always exhibits some part of randomness, 
either due to internal reasons, specific for nonlinear dynamical systems, or caused by 
external stochastic noise. Therefore randomness is inavoidable to this or that extent 
and, for some time intervals, it can mask the existence of underlying tendencies and 
persisting trends. Nevertheless, such trends arc to be noticeable on average, whether 
this concerns the evolution of financial markets or biological species. Persisting features 
retained during this evolution imply that the complex system evolves self-similarly. It 
is worth emphasizing that the Law of Self-Similar Evolution does not mean and in 
no sense requires that there should be displayed some rigidly fixed properties of the 
considered system, but it rather tells that there must exist some trends in evolution. 
The evolutional self- similarity is not a static notion requiring the occurrence of station- 
ary fixed properties but it is a dynamic concept stating the persistence of underlying 
trends. 

In this section, the Law of Self-Similar Evolution has been described as a philosoph- 
ical category, in rather general and perhaps vague words. The corresponding mathe- 
matical realization will be presented in the following sections. The suggested approach 
is mathematically based on the Self-Similar Approximation Theory whose application 
to the analysis of asymptotic series is expounded in Sec. 2. The self-similar extrapola- 
tion of asymptotic series can be directly reformulated as a self-similar forecasting for 
time series, which is explained in Sec. 3. A method for evaluating the probabilities 
of self-similar patterns is developed in Sec. 4. These sections present, for the ffist 
time, the complete mathematical foundation of the self-similar approach to analyzing 
arbitrary time series. Specifications related to market time series are also discussed. 
Principal points of the approach are summarized in Sec. 5. 

2. Self-similarity in asymptotic series 

Assume that we are interested in finding a function f{x) of a real variable x. Let 
this function be defined by so complicated equations that we are able to extract from 
them only asymptotic expansions in the vicinity of some point Xq. Without the loss of 
generality, the expansion point may be taken as zero, Xq — 0. Let us have several such 
expansions. 



enumerated by the index A; = 0, 1, 2, . . .. The basic problem is what can be said about 
the value of f{x) at finite x if all we know are asymptotic expansions fk{x) in the vicinity 
of x — > 0? This problem is constantly met in physics and applied mathematics, where 
it is often called the problem of function reconstruction or the problem of summation 
of asymptotic series [12]. 



[10,11]. 





(1) 
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An original general approach, named the Self-Similar Approximation Theory, for 
reconstructing functions from a set of their approximate expressions has been devel- 
oped and successfully applied to various problems [13-23]. The name comes from the 
basic idea of the approach to present the passage between subsequent approximations 
as a self-similar transformation. More precisely, it was shown [15-19] that for an ap- 
proximation sequence {fk{x)} it is possible to construct a cascade, that is a dynamical 
system with discrete time, whose trajectory is bijcctive to the sequence {fk}, so that 
the sought function f{x) corresponds to a fixed point of the cascade. If we treat the 
given complicated equations, together with a calculational algorithm, as a complex 
system generating a sequence {fk}, then the latter is nothing but a prototype of a 
time series, the approximation number k playing the role of time. Hence the self- 
similar approximation theory is the mathematical realization of the law of self-similar 
evolution. 

Asymptotic expansions are usually presented as power-law series 

fkix) = E «n ' (2) 

n=0 

where are arbitrary real numbers arranged in ascending order, 

< (n = 0, 1, . . . , A;) . (3) 

For the purpose of self-similar analysis, the presentation of the series (2) has to satisfy 
several general properties: 

(i) Analysis should not depend on the choice of units for the variable. Hence the 
latter is to be taken in a dimensionless form. 

(ii) Asymptotic expansions are to be reduced to a scale-invariant form. To this end, 
it is always possible to factor out the term 

fo{x) = ao x^° (ao 0) (4) 

and to introduce the scale-invariant function 

fo{x) 

where x is assumed to be dimensionless. Evidently, function (5) does not depend on 
the change of scales for fk- 

(iii) If the variable x pertains to a finite interval, the latter is to be normalized to 
the unitary interval. So that everywhere in what follows, it is assumed that x e [0, 1]. 

Note that if the variable pertains to an infinite interval, one may pursue different 

ways depending on what additional information on the behaviour of f{x) at a; — cx) is 
available. If there is no such information, it is possible to transform the infinite interval 
to the unitary one by means of the change of variables x' = x/{l-\-x), x = x'/ {1 — x'), 
so that x' e [0, 1]. 

For the scale- invariant function (5), we have 

k 

n=0 
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where x e [0, 1] and 

bn = — , /3„ = an - ao > /5o = . (7) 

Clearly, (fio{x) — 1. Also, all powers /3„ are positive, even if some q;„ are negative, 
which follows from the ascending order (3). Recall that, by construction, the scries 
(6) are assumed to be asymptotic, having sense only for a; — 0, while for any finite x 
sequence {ipk{x)}'^Q diverges. It is possible to say that the latter sequence converges 
just at one point x — 0. 

To proceed further, we have to transform divergent series to a form that would have 
sense for finite x G [0,1). This can be done with the help of control functions [13], 
which can be introduced in different ways [13-23]. Dealing with asymptotic series, it 
is convenient to invoke the multiplicative power-law transformation [20-22] defined as 

^k{x,s) = x' ipk{x) , (8) 

with the inverse transformation 

ipk{x) = X-' ^k{x,s) . (9) 

Since power laws are common in describing fractal objects, equation (8) may be called 
the fractal transformation. The transform (8), according to equation (6), is 

k 

^k{x,s) = Y^ hnX'^f"- . (10) 
n=0 

Here s = s{x) is a control function, whose role is to make the series (10) meaningful for 
finite X. These series can be considered as an expansion in powers of the new variable 
x^. As is evident, such series are asymptotic with respect to x^ 0. The latter limit 
can be achieved if, instead of forcing x to zero, we keep < 1 and setting s oo. 
Thus, the series (10) can be treated as asymptotic with respect to s ^ cxo for all |a;| < 1. 
Now we may say that the sequence {$^(2;, s)}^q converges for all \x\ < 1, provided 
that s — > 00. In this way, we come to the natural choice of the control function s — > 00. 
Recall that we are considering the case when no additional constraints are imposed on 
the behaviour of the sought function and all we know are its asymptotic expansions 
(1). In the intermediate expressions, the value of s is assumed to be asymptotically 
large, and the actual limit s — > cxo is to be taken after the inverse transformation (9). 

Since our consideration here concerns functions, we need to define the property of 
functional self-similarity, which should not be confused with geometric self-similarity 
describing fractals. The notion of geometric self-similarity [24] is connected with the 
scaling of a variable, which is only a particular kind of the more general notion of the 
functional or group self- similarity [15-19]. To correctly define the latter, we need to 
introduce some notation. We define the expansion function x{(p, s) by the equation 

%{x,s)=ip, x = x{ip,s). (11) 

With the form (10), this gives x{ip, s) — cp^^^ and cp — x^. Introduce the mapping 

yk{(p,s) = ^k{x{(p,s),s) . (12) 
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Let ^k{x, s) be real for all A; = 0, 1, 2, . . . and all x G [0, 1]. From ^k{x, s) G R it follows 
that yki^p, s) G R for all A; = 0, 1, 2, . . . and (/? G R. Therefore the mapping (12) is an 
endomorphism on R. It is this endomorphism that allows us to formulate the property 
of group self-similarity we need. 

Our aim is to present the change of the endomorphism (12), when varying the 
approximation number k, as the evolution of yk with respect to the discrete time k. 
From the point of view of group theory, self-similarity is nothing but a semigroup 
property. The latter, for the evolution of yk with respect to k, reads yk+p = yk ■ yp- As 
follows from the definition (11), the unit element is yo, since yo{^, s) = ip. The family of 
endomorphisms, {yk\ k — 0,1,2, . . .}, with the semigroup property forms a dynamical 
system in discrete time, called the cascade. The semigroup property, in terms of the 
notation (12), takes the form 

Vk+pi^P, s) = VkiVpi^P, s),s) . (13) 

Since in the accepted interpretation, we treat the sequence {yk{f,s)} as a trajectory 
resulting from the evolution of the cascade {yk\ k — 0,1,2, . . .}, the relation (13) may 
be called the evolutional self- similarity. As far as the corresponding semigroup property 
is natural for dynamical systems, the equation (13) may also be termed the dynamic 
self- similarity. This equation (13) is a necessary condition for the fastest-convergence 
criterion [16,17]; the cascade fixed point representing the sought function. 

Following the general theory [14-19], the cascade {yk\ k — 0, 1, 2, . . .} can be em- 
bedded into a flow {y-rl t > 0}, which is a dynamical system in continuous time. For 
the latter, one may write the Lie equation which is a differential equation of motion. 
The flow velocity field is defined, by means of the Euler discretization, as the cascade 
velocity 

Vn{(p,s) = yn{(p,s) -yn-i{(p,s) {u ^ 1,2, . . . , k) . (14) 
From equations (10) - (12), it follows that 

k 

yki^,s)^Y. bnip'^"-/', (15) 

n=0 

which results in 

r;„(^,s) = 6„/+^"/^ (16) 

The differential equation of motion can be integrated. The integration over the effective 
time goes from t — ntOT — n-\-Tn, with t„ being the effective time required for reaching 
a fixed point after the n-th step. In this way, 

fn+T„ 

/ dt = Tn. (17) 

Jn 

Therefore, the evolution integral [15-19] acquires the form 

dip , , 

This, with the cascade velocity (16), yields 
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where = ^n{x, s) = yn{x^, s) and = ^n{x, s). In particular, 

Returning back to the function ^pl{x) by means of the inverse transformation (9), 
we have to take the hmit s — > oo. For example, 

(pl{x) = hm x~^ ^lix,s) , 

which results in the first-order self-similar approximant 

(p*i(x) — exp (biTi x^^) . 

To obtain higher approximations, we can accomplish such a renormalization procedure 
2k times for each ^l{x, s), which was called [21,22] self-similar bootstrap. However the 
same result can be reached twice faster, by accomplishing k renormalizations, in the 
following way. We may present the function (6) as 

cpJx) = l + b,x^' + ^ x^'-^' f 1 + x^^'-^' (1 + • • •) 

The latter can be written in the form (pk{x) = 1 + Xi, in which Xi is expressed through 
X2, and X2 through 2:3, and so on according to the rule 

T — ^/3n-/3n-l n I ^ \ 

where n = 1,2, ... ,k. Considering each Xn as a small parameter, we need to accom- 
plish k times the first-order renormalization procedure described above. As a result, 
introducing the notation 

c„ = -^T„, Un = an-an-i (n = 1, 2, . . . , /c) , (19) 

we come to the fc-order self-similar exponential approximant 

ipl{x) = exp {cix"^ exp {c2x''^ ... exp {ckx'''^)) . . .) . (20) 

This, for short, can also be named the A;-order superexponential. 

It is worth emphasizing that although the form (20) reminds the Euler nested 
exponentials [25,26], it is principally different from the latter. First of all, the Euler 
superexponentials are defined only for integer powers «„. Second, when one tries to 
sum power series by means of such continued exponentials, one fits the coefficients 
in the latter so that to reproduce those in the power series, as a result of which the 
constructed superexponentials have the same radius of convergence as the related power 
series on the real axis [27,28] . 

In our case, the self-similar exponential (20) contains the coefficients c„ given by 
equation (19), where the effective renormalization time t„ is yet undefined. The latter 
plays the role of a control function that is to be determined from additional conditions. 
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One way could be to find r„ from tlie fixed-point conditions [22] having the form of 
the minimal-difference criterion [13]. Such fixed-point equations not always possess 
solutions. Here wc suggest another, and more general way of defining the control 
functions r„. This method is commonly employed for defining control functions in the 
optimal control theory. The idea is to construct a cost functional whose minimization 
yields the control functions of interest. 

The quantity r„ appears in the evolution integral (18), where it has the meaning of 
an effective time required for reaching, after the n-th step, a fixed point representing 
the sought function. One, clearly, would like to reach the answer as fast as possible. 
The minimal number of renormalization steps is, obviously, one. Therefore one would 
like that the total effective time nr„ be also close to one. Hence, one should look for t„ 
being close to 1/n. At the same time, how fast one reaches the fixed point depends on 
the distance of the latter from the starting point. The distance that is passed during 
the time r„, with a velocity f„, can be evaluated as f„r„. Thus, we need to find a 
minimal time t„, being close to corresponding to the minimal distance VnTn- This 
suggests us to construct the time- distance cost functional 



2 V 



1\2 

Tn + (VnTn) 

nj 



(21) 



whose minimization with respect to the control function t„ yields 

1 



n{l + vl) ■ 



(22) 



The velocity v„ — Vn{x) is to be understood as the image, in the domain of x, of the 
cascade velocity (16), with taking account of the inverse transformation (9), which 
gives Vn{x) = x~'^Vn{x^ , s). Prom here, one has 

Vn{x) = Lpn{x) - (Pn-l{x) = 6„ x'^" , (23) 

which defines the control function r„ = Tn{x) according to the expression (22). 
In this way, we find the controllers 

rn{x) = 2 2( ^^ ' ^n(^) = - • (24) 

n[l + vf^[x)\ ao 

Combining equations (5) and (20), we obtain the /c-order self-similar exponential ap- 
proximant 

fl{x) = fo{x) cxp {cix"' exp {C2x''' . . . exp {ckx"")) . . .) , (25) 

in which c„ = c„(a;), and fo{x) is given by formula (4). The approximant (25) extrap- 
olates the asymptotic series (2), valid only for x — > 0, to the region of finite x e [0, 1). 
The value of f{x) at the point x — 1 can be defined as the limit from the left, as 
X — > 1 — 0. Therefore the superexponential (25) extrapolates the sought function f{x) 
from asymptotically small x — > to the whole unitary interval < x < 1. 
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It is worth recalling that in the asymptotic series (2) the powers were assumed 
to be arbitrary real numbers, with the sole requirement that they are arranged in 
the ascending order (3). Some, or even all, of these powers could be negative. If so, 
then the initial term (4), with ao < 0, has the power-law form that has been so much 
discussed in literature. The asymptotic existence of power laws is well known in critical 
phenomena. The relevance of power laws has repeatedly been claimed to describe many 
natural phenomena, ranging from earthquakes [24,29,30] to different economic and 
financial distributions [31,32]. Since such power laws are practically always asymptotic, 
their more general form should include corrections leading to the power-law series (2). 
Extrapolating these asymptotic series in the described self-similar way, one should come 
to the self-similar exponentials. The first-order approximation then results in a kind 
of a stretched exponential. The stretched exponential distributions describe many 
phenomena in nature and economy either not worse or even better than power-law 
distribution functions [33] . More generally, the extrapolation of power laws should lead 
to the self-similar nested exponentials (25), whose structure evidently demonstrates the 
existence of many scales. 

3. Self- similarity in time series 

The technique of the self-similar extrapolation for asymptotic series can be refor- 
mulated as the method of forecasting for time series. It is necessary to call attention 
to the existence of several principal points in this reformulation. Overlooking these 
points would essentially restrict the applicability of the method. The correct general 
way of self-similar forecasting for time series is advanced below. 

First of all, the same basic requirements that were imposed on asymptotic series 
are compulsory for time series: 

(i) The measured quantity is to be presented in a scale- invariant form. This can be 
easily done by normalizing the given data to the value /o of the measured quantity at 
the initial time, which is analogous to introducing the scale- invariant function (5). 

(ii) The time variable has to be normalized to a dimensionless form, such that 
the prediction time would pertain to the unitary interval [0,1]. This normalization 
eliminates the ambiguity in defining the power of the power- law transformation (8), 
requiring that s 0, which results in self-similar superexponentials (20). 

In addition to these requirements, common both for asymptotic as well as for time 
series, there arises a question specific for time series: how the latter should be pre- 
sented, as a backward or forward recursion? The answer to this question follows from 
the comparison of asymptotic series with time series. In order that the extrapolation 
of asymptotic scries could be directly extended to forecasting for time series, the pre- 
diction horizon is to be a unitary interval and the available information from the past 
has to provide approximate forecasts for asymptotically small time t ^ +0. This re- 
quirement can be strictly accomplished only for backward recursion. Hence, we have 
one more restriction, specific for time series: 

(iii) For the correct usage of the self-similar analysis, time series are to be arranged 
as a backward recursion. This means that, if we are given a set {/„} of data corre- 
sponding to a quantity of interest, measured at the times tn, where n — 0,1,2, ... ,k, 
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the moments of time are to be ordered so that 

tn+i<tn {n = 0,l,2,...,k) , (26) 

with the initial time to — 0. The past-history data base of k-th order is the set 

Dfc = {fk, fk-i, ■■■,fo\tk< tk-i < . . . < 0} . (27) 

To start the procedure of self-similar forecasting for future times t G [0, 1], we need 
to possess a sequence of approximations valid for asymptotically small t — > -|-0. The role 
of such asymptotic forecasts can naturally be played by functions fk{t) interpolating 
the given data base (27) for the past time horizon tk < t < 0, so that 

fk{tn) = fn {n = 0,l,...,k) . (28) 

This interpolation can be uniquely defined by the Lagrange interpolation formula [34] 
presenting the interpolation function 

k 

fk{t) = Y. fn Init) {k > 1) (29) 

n=0 

as a series over the Lagrange polynomials 

Init)^ n 1^ in<k). (30) 

Because of the property Initm) = ^mn of the Lagrange polynomials, condition (28) 
is automatically satisfied. Note that the interpolation form (29) is evidently scale 
invariant with respect to time and the ratio fk{t)/ fo is scale invariant with respect 
to units of /fe. Also, one may remark that the time moments t„ are not necessarily 
equidistant, but can be chosen arbitrarily. 

The interpolation formula (29) can be rewritten as the algebraic polynomial 

k 

fk{t) = E «n (ao = /o) , (31) 

n=0 

with the coefficients a„ = a„fc immediately following from equations (29) and (30). The 
values of these coefficients depend, of course, on the given data base (27), but for the 
simplicity of notation, we shall omit in what follows the additional index k. Consider 
the subsequence {fn{t)}n=o ^f the terms 

Ut) = X: «m , (32) 

m=0 

which, for a given k, tends to the polynomial (31) as n — > /c and for which 

/„(0)=ao = /o (n = 0,l,...,A;). (33) 

Then the subsequence {fn{t)}n=o '^^^ ^c treated as a sequence of approximations 
asymptotically vahd for t — > -|-0. Employing the self-similar extrapolation, described 
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in the previous section, we may construct a forecast for the time interval < t < 1. 
The result is given by formulas (24) and (25), where we need to change the variable 
X by the time t and to put Q!„ = n, = 1. Instead of equations (24), we have the 
controllers ^ 

Cn{t) = — ^ Tnit) (n = 1, 2, . . . , A;) 

= -iv^rm ' = T ■ (34) 

And the self-similar exponential (25) takes the form 

fkit) = fo exp {cit exp {c2t . . . exp (ckt)) . . .) , (35) 

in which Cn = Cn{t). Expression (35) is the self- similar forecast for the future time 
interval <t < 1, predicted on the grounds of the past-history data base (27). 

4. Ensemble of possible scenarios 

For each given data base (27), one may construct the self-similar forecast (35). But 
one can take several different data bases by varying either the data-base order k, or by 

changing the time intervals between the moments of time t„, where n = 0, 1, . . . , A;, or 
by varying both the data-base order as well as the data-base scale. So that, in general, 
one can consider an ensemble of different data bases. Each of the latter has to be 
labelled by two indices, 

D.(j) = {fi'\ f!f}„ . . . , /ol ^ ^ < ^2, < . . . < 0} , (36) 

one index, k, defining the data base order, and another, j, specifying the chosen time 
scale of the past. For each data base (36), one obtains a forecast f^{j, t) according to the 
rule (35), but with different values of c„ = Cnk- Hence there exists an ensemble {/^ (j, t)} 
of possible forecasts, or admissible scenarios. Which of these possible forecasts should 
one trust? 

Life teaches us that in the majority of cases nothing can be trusted for hundred 
percent. But, when there can happen several different events, they can be classified by 
estimating their probabilities. Hence we need to define a probability measure on the 
ensemble of scenarios {/^(j, t)}. 

This problem is analogous to the problem of pattern selection occurring for non- 
linear differential equations in partial derivatives. Such equations sometimes possess 
a set of solutions corresponding to different spatio-temporal structures, or patterns 
[35]. A general approach for treating the problem of pattern selection has been sug- 
gested [36]. This approach can be directly applied for weighting possible scenarios 
from the given ensemble of self-similar forecasts. For this purpose, we may consider 
the passage from /^(j, i) to /fc+i(j, t) as the motion with respect to k. Then the map 
{fk{j, t) I /c = 1, 2 . . .} is to be treated as the image of a dynamical system with discrete 
time k. The probability of a scenario fk{j,t) can be defined as 

Pk{j, t) = ^ exp {-ASkU, t)} , (37) 

^k(t) 
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where Zk{t) is a normalization factor, being the sum 



Zfe(t)=^exp{-A5fe(j,t)} 



over the pattern indices j, and the entropy variation 

AS,{j,t) = S,{j,t) - S,{j,t) 



(38) 



shows the change of entropy with respect to the effective time k. The entropy of a 
dynamical system may be defined, by analogy with statistical systems, as the logarithm 
of an elementary phase volume [36], the latter, in our case, being \Sf^{j,t)\. Thus the 
dynamical entropy is 

S,{j,t)^\n\dflij,t)\. (39) 
Then the entropy variation (38) becomes 



ASk{j,t) = ln 



(40) 



With the notation for the mapping multiplier 



dmj,t)/dt 



the entropy variation (40) reduces to 

A5'fe(j,t) = In \mk{j,t)\ . 
It is convenient to introduce the average multiplier mk{t) by the relation 

1 ^ 1 



|"^fc(^)| 



\mk{3,t) 



(41) 



(42) 



(43) 



Using equations (42) and (43), for the scenario probability (37), we have 

mk{t) 



mk(j,t) 



(44) 



Note that the scenario probability (44) is normalized with respect to the summation 
over the pattern indices j corresponding to different data-base scales. In a particular 
case of just one fixed scale, one could accomplish the normalization with respect to the 
summation over k, which would define the probability weights for a restricted data- 
base ensemble, as was postulated in reference [37]. The derived scenario probability 
(44) concerns the general case of an arbitrary ensemble {Dfc(j)} of the data bases (36). 
Being general, the approach of the present paper makes it possible to answer several 
principal questions. 

One important question concerns the choice of the data-base order. More in detail, 
the problem is as follows. For a fixed time scale, labelled by j, we may analyze different 
data bases Dfe(j), with varying /c = 1, 2, 3, Then, how many terms fjf^ should we 
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take? That is, when should we stop increasing kl The answer is straightforward: 
The data-base order k has to be increased till we reach numerical convergence. More 
precisely, this means the following. Let us be satisfied by the results of an error s. 
Then we need to increase k up to the saturation number Nj = Nj{s) such that 

\fUnU,t) - rkU,t)\ < ^ (45) 

for k > Nj, all n > 0, and t e [0, 1]. For this k — Nj, we have the saturated data base 
T>{j) = T>N.{j) and the related saturated forecast 

nj,t)^r^^u,t) (46) 

characterized by the saturated scenario probability 

p{j,t)=pN,{j,t) . (47) 

Thus, varying the data-base scales, labelled by j, we get the saturated ensemble {f*{j, t)} 
of scenarios (46), with the probability measure (47). 

The most probable scenario from the ensemble {f*{j,t)} is the forecast f*{jo,t) 
having the largest probability, such that 

maxp(j, t) = p{jo, t) . (48) 

The latter, because of the form (44), is equivalent to the condition of the minimum for 
the absolute value of the multiplier m{j,t) = mp^.{j,t), so that 

min |m(j, t)\ = |m(jo, ^)| ■ (49) 

3 

Having the probabihty measure (47), it is possible to define the average forecast, or 
the expected forecast 

<f{t)>^j:p{j,t)f*{j,t). (50) 

3 

The dispersion 

a'{t) ^ < f\t) > - < fit) , 

in the case of a market, describes the market volatility. The latter can also be charac- 
terized by the variance cocfiicicnt K,(t) = o'(t)/ < f{t) >. 

The last question to be answered in order to have a completely self-consistent 
theory is how to choose the data-base time scale. In general, the moments of time t^^^ 
should not be compulsory equidistant. However, for practical purpose, it looks more 
convenient to take them as such, defining the time step as = t^p — t^^^i. It seems 
natural to start with the time scale Aq = 1, which equals the prediction horizon. Then, 
one may decrease as well as increase the time step, for instance, according to the rule 

= 2~^, ^2i+i = 2-^ , with J = 0, 1, 2, One has to stop decreasing and increasing 

the data-base time scale at such j — jmax, when numerical convergence is reached. 
This implies that, for a given error e, one gets the inequality 

\f*(j+m,t)-f*{j,t)\<s 
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for j > jmax, all m > 0, and t e [0, 1]. 
5. Summary of main ideas 

In this paper, a novel approach to analyzing time series has been presented. There 
are several principal points distinguishing this approach from the standard one. The 
aim of this paper has been to clearly describe these principal points forming the body 
of a self-consistent theory. Not yet all parts of this theory have been exploited in full 
for practical applications; the calculational work is in process. But some simphfied 
versions of the approach have been illustrated by a number of examples for market 
time series [37,38]. In this concluding section, I would like to emphasize again the 
main ideas the approach is based on. 

To better stress the principal difference of the present approach from the standard 
way of analyzing time series, let us recall the basic idea of the latter: For a given set 
of data {fn\ n — 0,1,2, . . . , k} one tries to invent a relation ft — /(/o, fi, ■ ■ ■ , fk, 
connecting the value ft at the moment of time t with the past data. This relation 
can be in the form of an explicit function or in the form of a difference or differential 
equation, including a stochastic term ^ modelling noise [1-5] . There are, to my mind, 
two principal deficiencies of such an approach. First, I think that no explicit equations, 
indifferently to how elaborated they are, can grasp all peculiarities of a realistic complex 
system, such as a market. Second, any given relations reflect only the past history, 
providing an interpolation for the learning historical period, while for predicting future 
one needs an extrapolation. 

The self-similar approach is based not on attempts to invent a relation between the 
historical points but it tries to discover dynamic trends resulting in these points. For 
this purpose, instead of studying relations between points, it is necessary to analyze 
relations between sequences. Interpolative formulas are used here only as a starting step. 
The comparison of different interpolative expressions makes it possible to produce an 
extrapolation, that is, forecasting. 

The basic philosophical idea of the self-similar analysis of complex systems is the 
Law of Self- Similar Evolution, formulated in Sec. 1. And the mathematical foundation 
is provided by the self-similar extrapolation of asymptotic series, described in Sec. 2. 
This extrapolation can be reformulated as forecasting for time series, as is done in 
Sec. 3. Since any forecasting can only be probabihstic, the way of defining the related 
probability measure is explained in Sec. 4. Some details of the self-similar analysis can 
be changed. For example, the number of points in a given data base can be reduced 
by replacing several neighbouring values by either arithmetic averages or by fitting 
the values on a large time interval with the help of simple splines [34]. However, in 
the process of these kinds of averaging, some information on the considered time series 
will be inavoidably lost. 

Another change could concern the definition of the effective control time t„. The 
latter is to be defined from the minimization of a cost functional. As an example, 
the time-distance cost functional (21) was considered. But, in general, one could 
opt for another cost functional, depending on the available information and imposed 
constraints. 

Also, it would be possible to deal not with the initially given data base but with 
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some transforms of it. For instance, one could consider the set {ln/„} instead of {/n}- 
Or one could keep in mind a more elaborated transform, like the wavelet transforms, 
often employed for analyzing time scries [39]. 

It is worth stressing once more that the notion of self-similarity exploited through- 
out the paper is understood here as the group self-similarity, which is a more general 
notion than the trivial geometric self-similarity one usually talks about in connection 
with fractals. In the latter case one assumes the existence of the scaling relation 
f{\x) = A"/(,x) for the considered function f{x). Such a relation, with a given bound- 
ary, or initial, condition f{xo) = /o, immediately results in the power-law function 
f{x) — fo{x/xo)°'. One could consider a slightly more complicated scaling relation 
as f{Xx) — u{X)f{x), with a known function u{X). However again, with the given 
boundary condition, this immediately gives the answer f{x) = fQ/u{xQ/x). All such 
scaling relations produce the considered function in an explicit form. Whereas the 
group self-similarity (13) provides an equation that is yet to be solved. 

An important point is that the group self-similarity employed here, and which is the 
basis of the Self-Similar Approximation Theory [13-23], has to do not with a scaling of 
a variable but with the motion with respect to the effective time whose role is played 
by the approximation number. Here it is the motion on the manifold of approximants. 

Finally, if we would like, digressing from mathematical foundations, to briefly con- 
clude why the idea of the group self-similarity does work for extrapolating asymp- 
totic series and forecasting time series, then we should return back to the Law of 
Self-Similar Evolution telling that all complex systems develop self-similarly, preserv- 
ing their generic features in the course of their evolution. Such features may be not 
noticeable from the first glance but rather hidden somewhere in genes. The group 
self-similarity is a kind of genetic self-similarity. Remember also that the Lord created 
the man in a self-similar way [40]. 
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